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We study a particle on a ring in presence of a dissipative Caldeira-Leggett environment and 
derive its response to a DC field. We show how this non-equilibrium response is related to a flux 
averaged equilibrium response. We find, through a 2-loop renormalization group analysis, that a 
large dissipation parameter n flows to a fixed point rj R = h/2n. We also reexamine the mapping of 
this problem to that of the Coulomb box and show that the relaxation resistance, of recent interest, 
is quantized for large n. For finite r\ > r\ R we find that a certain average of the relaxation resistance 
is quantized. We propose a Coulomb box experiment to measure a quantized noise. 
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I. INTRODUCTION 

Two of the most important mesoscopic structures are rings, for the study of persistent currents, and quantum dots or 
boxes, for the study of charge quantization. Of particular recent interest is the quantization of the relaxation resistance 
R q , defined via an AC capacitance of a single electron box (SEB). Following the prediction of Biittiker, Thomas and 
PretreS that R q = h/2e 2 for a single mode resistor, a quantum mesoscopic RC circuit has been implemented in a 
two-dimensional electron gas ^ and R q = h/2e 2 has been measured. The theory has been recently extended to include 
Coulomb blockade effects^' showing that R q = h/2e 2 is valid for small dots and crosses over to R q = h/e 2 for large 
dots. 

In parallel, recent data has observed Aharonov-Bohm oscillations from single electron states in semiconducting 
ringfP. Further theoretical works have considered the effects of dissipative environments on a single particle in a 
ringS, in particular studying the renormalization of the mass M* and its possible relation to dephasing^"^. A related 
case of a ring coupled by tunneling to an electron lead has also been studied^. 

It is rather remarkable that the ring and box problems are related via the Ambegaokar, Eckern, and Schon (AES) 
mapping^ where the ring experiences a Caldeira-Leggett (CL^l environment. While the exact mapping assumes 
weak tunneling into the box with many channels, it has been extensively used to describe various tunnel junctions^, 
the Coulomb blockade phenomena in SEB and in the single electron transistor (SET^P^H. 

The ring problem is defined by a particle confined to a ring, coupled to a dissipative environment of the Caldeira- 
Legget type, and in presence of a field E, generated by a time dependent flux through the ring. This scenario is 
schematically illustrated in Fig.l. In the present work we address the ring problem by the real time Keldysh method 
and study it using a 2-loop expansion and renormalization group (RG) reasoning. We find that perturbation theory 
identifies an unexpected new small parameter sin(^) where 77 is the dissipation parameter on the ring, or the lead-dot 

coupling in the SEB. We infer that a large 77 flows to a fixed point r/ R with h/2r] R — tt. Whi le the thermodynamics of 
the ring type problem has been much studied, including extensive Monte Carlo studies^^ of M* , no sign of a finite 
coupling fixed point has been detected. Our method evaluates the response to a strictly DC electric field E, equivalent 
to a magnetic flux through the ring that increases linearly with time, hence a non-equilibrium response. We claim 
that thermodynamic quantities like M*, that are flux sensitive, decouple from the response to E, a response that 




FIG. 1: Artist view of a particle on a ring, coupled to an environment, with a field E — cj> x due to a time dependent flux 
through the ring. The particle polarizes the environment which in turn modifies the motion of the particle at later times, i.e. 
an effective non-local interaction. 
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averages over flux values. This general relation between non-equilibrium and equilibrium responses is given by Eq. 
( 39 1 below. This relation has been noticed for a model with p arti cle tunneling between a ring and an environment^. 



In terms of the SEB, our results extend the previous analysis^' to the case of many channels N c , an experimentally 
realizable scenario^. We note that for N c > 1 the relaxation resistance for noninteracting electrons^ becomes 
h/(2N c e 2 ). We find that for strong coupling, 77/?! > 1 the relaxation resistance is quantized to e 2 /h up to an 
exponentially small correction ~ e - '" ? ' ft . p or finite 77, but still r\ > -q R we find that a certain average of the relaxation 



resistance is quantized (see Eq. 82 1 



The present work considerably expands our previous letter 2 ^. In section II we present the ring and box models, with 
some exact general properties. In section III we present RG and numerical solutions for the semiclassical case, while 
section IV presents the perturbation and RG analysis of the full quantum case. The discussion section V summarizes 
our results, presents a topological interpretation of our fixed point and details a proposed Coulomb box experiment 
to detect our predicted quantized noise. The Appendices give details of the ring-box mapping and of the various 
perturbation expansions. 

II. THE MODEL AND GENERAL PROPERTIES 
A. Semiclassical model 

We derive first a Langevin equation for a particle on a ring. Consider the standard Langevin equation for a particle 
with coordinate Xt in one-dimension of the form 

R- t }x t > = it (i) 

where £ t is a Gaussian random force from an environment, where the average on the environment degrees of freedom 
is 

(66') = B ttt , (2) 

This relation defines a linear response for either x u — R u £ u or £(cj) = i? -1 (w)ir(cj), after Fourier transforms, e.g. 
R(uj) is the Fourier transform of R t — Rt,o- Hence the fluctuation dissipation theorem (FDT) at temperature T can 
be applied either way, leading to 



K x (u) = ?jcoth(i ( 0/iu;)Im[.R 1 

-1 
R 



B u = ftcoth(|/3fio;)Im— (3) 



where K x (uj) is the Fourier transform of K x (t) = ^(xtXt+ T + Xt+T%t)- The simplest choice corresponds to a particle 
with mass m and a friction coefficient 77, so that at temperature T = 

mx t + r]x t = i t 

Ro(u) = ^— R (t) = -[1 - e~^ m }e(t) 

= h V \u;\ Bt = ^(t^0) (4) 

where Ro(t — t') is the response in this case. While the mass provides a high frequency cutoff which we denote 
uj c = -q/m, the singularity of B(t) at t = Q implies the need for an additional cutoff. This additional cutoff is a 
convenience and will be used below in the simulations as well as in the RG derivation. A method for deriving general 
response functions is based on Kramers Kronig relations^. In the notation of Eq. (2.7) of Ref. [3D] we choose 
Rc/j(u;) = 77/(1 + uj 2 Tq) so that the response function i? _1 (i — t'), after Fourier, is 

K, — —muj (5) 

which has the remarkable and necessary property that both and i?" 1 have no poles in the upper half plane. The 
FDT at T = gives 



(6) 
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so that 1/ro provides a cutoff on the environment frequencies. Hence for 4to < m/rj, {8 — > +0) 

TO 77T0 AiX XiX 

\x = \[l+x], A 2 = §[l-4 x =Jl-^T° (7) 



while for 4r > 771/77 with x = y ^f- - 1 

i?t = 0(t)-{e~ 5t - [i^- sin(^/2r ) + cos(rf/2T )]e-*/ 2ro } (8) 
77 2x 

Consider now the two-dimensional system and its projection on a ring, i.e. x t = (cos 8 t , sin 9 t ) so that # t is the 
angular position of the particle and the radius is chosen as unity. In cartesian coordinates we define random forces 
in the x,y directions so that R^_ t t cosdf — ~Cti sm = &• The ring potential confines the motion to the 

azimuthal part, so that only the tangent force — £f cos 6^ + £ f b sm9 t is allowed, hence 

- sin 6(t)R-\, cos B t , + cos 9 t Rjl v sin 9 t , = cos t + $ sin 6 t + E (9) 

where Q, $ are independent and each having the correlations of Eq. An external tangent electric field E has 
been added corresponding to a flux through the ring that is increasing linearly with time <p x = Et. With Ro(t — t') 
given by Eq. Q the differential form Rq 1 ^) — mrdf + rjrdt, can be used leading to 

md\ + T)6 t = Q cos 9 t + $ sin 6 t + E (10) 

Comparing the terms on the left identifies a cutoff frequency uj c — rj/m. At u) > lj c the mass term dominates while at 
lj < u> c the environment dominates, leading to renormalizations. The nonlinear Langevin's equation (10), including 
an average on the random forces, is equivalent to a partition function 

Z = J V[6,€\ 8 (m9 t + rj9 t - cos 8 t - $ sin 6 t -E) e" IJ\tZ\*+K\ 3 V 2B » 

(11) 

Introducing the 'quantum' field 6 by 8(X t ) = J V[9]e ietXt , and averaging over the noise field £ x , £ y results in the semi 
classical partition function Z = J V[9, 9]e~ s ^ e ' s ^ where S[9, 9] — So + Si n t is given by 

S = i [ § t (Rt, t .)-% -*e [ e t ,=i[ R^e u 0- u - 1E 

Jt,t' Jt' Jul 

Sint = \ I o t B t ,tA cos{e t -e v ). (12) 

h,t' 



2 



This has the form of a Keldysh action, with 9, 9 being the classical and quantum fields, respectively. We will see 
below that this action is the semiclassical H — > limit of the full quantum system. 



B. Quantum model 

We proceed to define the full quantum problem. The one-dimensional Langevin systenf^^ has the Keldysh 
partition Z = J r Dx t T>x t e~ SK where 



Sk =i / x t R t },x t , + | / x t B, 
Jt,f Jt.t' 



t.t'Xf (13) 



and Xt, Xt are the quantum and classical fields, respectively, 



x t = \{xt + x t ) x t = -(xf-x t ) (14) 
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and xf 1 are on the upper and lower Keldysh contour, respectively. On a ring, we use a 2-dimensional vector notation 



(15) 



Denning 



and using trigonometric identities we obtain the quantum action 

S K=ilj Rtj> sin(^t) cos&) sin(0 t , - t ) + 4 / smfy t ) sin&) cos(0 t , - 9 t ) 

J t.t ™ J t.t " 



(17) 



We note that the path integral involves continuous 9 t trajectories that can involve n rotations around the ring. 
Consider the time evolution from an initial wavefunction ipifio, to) at time to to a final state tp(9t, t), where both initial 
and final angles are compact, < 9 , 9 t < 2ir, 



il>(8 t ,t)= d6 y2 Vee- s ^iP(6 ,t ) (18) 

Jo n Je 

The sum on the integers n expresses that the probability to arrive at a given 9 t is a sum of probabilities, each with 
n rotations. The path integral can therefore be written in terms of a decompactified variable 9 t =_0t + 27m, i.e. 

Sn le t+2 ^ n 'D® ~~ ^ le* where now — oo < 9 t < oo. This shift does not affect the periodic forms in (17), however it 



does affect an external electric field E. Consider a time dependent flux <p x (t) — Et that contributes to the action a 
term J f <f> x (t)6 t dt — —E J t f 9 t dt + (f> x (t i )9t i —</> x (tf)8t f . The partial integration is allowed only for the decompactified 
variable 9 t , i.e. the work done by E is finite for each 2tt rotation. The boundary terms are neglected, e.g. one can 
choose 4> x {ti) — 4>x{tf) — where ti,t/ — > — oo are boundary times on a Keldysh contour; the field E is turned on 
slowly away from these times. 

In the following we will consider a perturbative scheme with a field E and a bare velocity v — E/rj and 9t is 
decomposed to 9t = 59 1 + vt; (the true velocity is defined below as v R (E) — (9 t )). The velocity v provides a low 
frequency cutoff eliminating divergence of the perturbative expansion and eventually allows for RG treatment. It will 
be convenient to use the two-cutoff response Eq. |H| with R^ 1 = —mux 1 + SR^ 1 , where SR^ 1 = jE^^j hence 

SR-}, = d t > [ ^ B -Mt-f) = _?L dt J e -(t-t')/T 0e r t _ t >y = !L e -(t-t')/T 0e , t _ t ^ dt/ (19) 

J u l- IUJT T T 

The operator identity is satisfied for any function decaying faster then e'* '/ r ° at t' — > — oo. Note, 

9 t 5R-l,vt' =i— fe t [ e- (t - t,)/To dt' = ivri f 9 t (20) 



t,t , 

t,t' T Jt J-oo Jt 



The mass term with mw 2 — > 5(t — t')dtdt> produces m J. 9 t 9t — m L 9 t S9 t + mv J. 9 t ; the last term with mv = E/l*j c 
is neglected relative to the field term J. Et9t- The full action is then 

-S e 

f 9 t R~}9 t , - iE [§ t = i f 9 t R u }59 tl 
Jut 1 Jt Jt,v 

Bt,f sin(^0 t )sin(^0 t ,)cos(0 t , - 9 t ) 

6R-l,{sm(^9 t )cos(p t ,)sm(9 t , - 9 t ) - ^0 t v \ (21) 

The use of a single cutoff Eq. ^ with 

Ro 1 (t,t') = 6(t-t')[md t dt'+vdt>) (22) 



S K 


— So + S. 


Sq 


= i f § t 




Jt.v 






Sint 








So 


= it 
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leads to a simpler action. It corresponds to tq — > 0, hence SR t \, — \ r]S(t — t')df 



-Ro 1 ^, t') sin(^ t ) cos(^) sin(0 t , - 9 t ) = 5{t - t')[m9 t 9 t + \ sm(h9 t )9 t - 



(23) 



where t is infinitesimal below t so that the retarded nature of R t t , is maintained. The action Sk = Sq + St n t + S c 
is then 

S a = i [ 9 t R^{t,t')59 t , = i [ [m§ t 50t + vOtSOt] = i [ [mOA + rfiA] -%E f 9 t 
Jt,t' Jt Jt Jt 

s 2 



, , Bt,v sin(-4) smfy t ,)cos(9 r - 9 t ) 
« Jt,f 2 2 



ii] 



sm(M t )9 t - -h9 t 9 t -\ 



r ^0. 



(24) 



Note that this action reduces to that to the semiclassical case Eq. (12) when h — > 



C. Renormalized friction 



The renormalized friction rj R (E) is defined by the renormalized response R R t , = i (e t §t>) n and its DC limit 



r] n (E) w-i-o 



(25) 



in analogy with the bare form Eq. (jij). We show now that the renormalized i] R (E) is also the local slope of 
where v R is the E dependent renormalized velocity 

Therefore 



V[9]9 t e- SK 



(26) 



dv R 

~dE 



I \ J 17(17,' 



pJL«, = / / H^c— (*"*') = lim = -j^™ 



(27) 



In particular we are interested in the limit rj R — r) R (E — > 0). 

We show now an alternative procedure for evaluating rj R . Consider the Keldysh partition Z = J T>[9]c~ Sk and shift 

9t —¥ 9t + at- The result must be at independent, and choosing the form (23) with tq — ¥ (the following identity is 
actually independent of cutoff choices) 



f) - 52 \ 
U - j— o 

oa t 



^(Sq + S int + S c ) s 



-i(rjv R -E-5E) 



./S(S int + S c ) 
6E — i 1 



(28) 



since -i(^r) — —m(9t) + r]{9) — E and v R is time independent, at least for long times. 



■ S9 t 



Taking an E derivative of Eq. (28) and using (27) we obtain 

1 1 Id 

r/ R (E) rj rj 1 dv 



SE 



(29) 



We have checked, up to 2nd order terms, that the results of (27) and (29) coincide. The use of (29) is technically 
easier. 
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D. Equilibrium correlations 

In this section we consider the equilibrium response to a change in flux and derive a relation with the nonequilibrium 
response to a field. 

Consider now the form of K(lj) as a response to a flux 4> x . Linear response to &li. r mg — +h9S<p x (t) is 

h{9) = - j ' K ttt ,8<t> x {t') (30) 
This corresponds also to the velocity correlation 

k t)V = +i9(t-t'){[9 t ,9 t >]) (31) 
We expect that the DC response is positive for small cj) x , hence define 

K{u) = -Kofo) + iuKifa) + 0{lu 2 ) (32) 



The response Ko{4> x ) is the persistent current, i.e. for a static flux one can integrate (30 1 

(9) = (' tf (C)#x (33) 



The periodicity of the persistent current implies K (4> x )d(f> x = 0. The curvature at tf> x = is a well studied object. 
For general <f> x it is defined by a Matsubara imaginary time correlation 

hW = ^ Jo Jo {WCdTdT ' = KM (34) 

where K(iuj n = 0) = +Kq (there is a sign difference in the standard Matsubara notation). An effective mass is defined 
by Ko(0) = h/M* so that M* = m without interactions, while for strong r\ ^> 1 coupling M* ~ e^ 17 is exponentially 
largePS. 

To appreciate the role of K\ consider FDT for the symmetrized correlation at small uj 

(%\ 2 ) sym = signo; ■ lmK u = \u\Ki (35) 
The diffusion involves the response (|#a;| 2 ) = Ki/\u\, hence for t — >• oo 

((9 t - 9o) 2 ) =Ki[ du 1 ' ^ = ^ ln( W ,t) (36) 

where uj x is a characteristic frequency where higher order terms in uj terms set in. 

Consider now the linear response to an electric field SH r ing = —E(t)6t and use the response (9t) — R^ t ,E(t') The 



definition (251 implies that the low u limit has the form Rfj = ^ ft . Since E = H4> x we expect ftw Rfj = K(to) 
However, t 



icrc is a difficulty with the latter relation, if taken literally. 



iujr) R 



l=?-K (<j> x ) + iuK 1 (ct> x ) (37) 



It is also not clear which <p x to use in this relation. To resolve this issue consider the K response with a constant 
electric field 

h(0 t ) = - I K t>v ■ Et 1 (38) 



Note first that an additional constant <j) x in \Et' + <p x can be eliminated by redefining the origin of the time if , 
hence the persistent current part should be eliminated. More precisely, define <j) x (t) = \Et\ the uj = component 
Koi'fix) = KoihEt) becomes a periodic function, i.e. an AC response with frequency u>e = ^E. For w — > this 

persistent current response averages to zero, i.e. J Q Ko(4> x )d(/) x = 0. The same reasoning applies to a cf> x average on 
Ki(<j) x ). Hence for the purpose of evaluating the DC response of (25) we need to average on the flux in (32), hence 

lim ii m = f Kri^d^ = \ . (39) 
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The order of limits in ^ signifies that t] R is essentially a non-equilibrium response. The equilibrium - nonequilibrium 
relation ( 39 1 has been noticed in solution of a Boltzmann relaxation equation for particles on a ring, allowing for 
particle tunneling into an environment^. 

The physical picture is that in a DC field the particle rotates around the ring and produces two types of currents. 
First is the persistent current that oscillates in time as <p x increases and is therefore time averaged to zero; this current 
is non-dissipative. Second, there is a genuine DC response from the iuKi term, which is dissipative. 



E. The Coulomb box 



Consider now the Coulomb box system, i.e. a finite region (a "dot") with charging energy E c coupled by tunneling 
to a single metallic lead. The Hamiltonian is 

H = e fc4,i a M + e ^ ati d a ,i + E C (N - N ) 2 + Y t k , a ,ia{ A d a ,i + h.c. (40) 

k a, i k,a,i 

where i = 1, N c are channel indices, d a ^ are dot electron operators with spectra e Q , a,k,i are lead electron operators 
with spectra e^, N = i i^a,i is the number operator on the dot, E c = e 2 /2C g is the charging energy with C g 
is the geometric (bare) capacitance, N is the gate voltage in units of 2E C . The channel index i is diagonal in the 
tunneling term, i.e. corresponds to transverse modes that are conserved in tunneling. 
Consider the density correlations 



K tfi = +i9(t-t'){[N t ,N t ,}) (41) 

The AES mapping to the ring problem is reproduced in Appendix A. In particular, Nq corresponds to — (f> x , 2E C to 
h 2 /m and the relation to the velocity correlation on the ring is 

h 2 K t , t , = -2E c h5(t - t') + 4E 2 K t , t , (42) 

Using the notatiorP K{uj) = H^§-(1 + iuiCoRq), where Co is the renormalized capacitance and R q is the relaxation 
resistance, we obtain 

hK(u) = -2E C + 4£ C 2 ^(I + iujCoRn) (43) 
Hence the mapping between the Coulomb box and the ring for the curvature is, using ( |34[ ) 

^n = ^o(«=2^(l-^) ^TT^TT-l-^ (44) 
M*((j) x ) C g M*(4> x ) Cg 

while for the dissipation, using p9] ) 

-4 = / KMd<l> x = % f C ^j^R q (N Q )dNo (45) 

V JO n JO 

We note that J* Co ^ No) dN Q = I due to the periodicity of F((f> x ). An extensive studjE^of M* (0) shows that it satisfies 
M*(0) > m and that for large rj (the bare interaction parameter) M*(0)/m ~ e m ' 3> 1. Hence 

^° = I - 0(6-") rj > 1 (46) 

and C — > C g for large r\. 

At this stage we can already propose an interesting experiment for the SEB. By analogy with E = hcj> x in the ring, 
we propose measuring the response to a gate voltage that is linear in time Nq ~ t. This leads to a DC current into 
the Coulomb box whose dissipation is the average in Eq. (45). This average is predicted to be quantized, at least for 
rj > r) R , as shown below. 



III. SEMICLASSICAL RG AND NUMERICS 



A. Perturbations and RG 



We study here the action (12) with a perturbation series for correlation functions. Consider first the correlation 



Cf,t — (Of ^t), which to 1st order is 

C { t \] = (9 t ,6 t (-S mt )) So = f B tuta convfa - t 2 )R tM R t 

Jti,t 2 



In Fourier space 



C^ = \R U \ 2 BZ 



(47) 



(48) 



where = \ (B U+VQ + B u _ Vo ). Since C t 9"2 t is divergent it is useful to evaluate C t ,t' — ([@t — 6f] 2 ), which to 1st order 
is, with t = t — f (t> 1/w c ), 



2h f lo S(^) T< -» 

C T = / £C|i? w | 2 (l-COSCJT) ' 



7T77 1 



(49) 



TTVT - < T 



For E — the angular position diffuses logarithmically, while for E ^ the long time fluctuation is linear in time. 
Consider next the response function to 2nd order in Si n ±, 



R", = i 



— Rt,f + R-tJi + Flf'}i — Rt,f + i (@t'Qt(—Sint + \Sf nt ) 



Sn 



(50) 



Note that the disconnected terms in the perturbation (S" nt ) Sq vanish for any order n, due to the normalization Z = 1. 
The first order response function is 



R { 3- = / b h m { K K ™s(e tl ~ e t2 )e v e t ; 
z Jti,t 2 x ' b ° 



i 



The result in frequency variable is (see Appendix B) 

i?« = Rl f R U1 [B V U1 - BZ_ U1 ] =Rlf R t B t cosv t (e iw * - 1) 

We note that for v = FDT is maintained, to this order, \v=o = lmR u h sign(w). 
The renormalized r] to first order is then 



(51) 



(52) 



1 



lim (^iui)Rj-^ = lim 



w->o {—ioj) 2 rf- 



RtB t cosvt (iojt) 



■log(l +u 2 c /v 2 ^ 



2rj' 



\ogv/u c 



■0{v) 



(53) 



Considering next the 2nd order in (50) we obtain (sec Appendix B) 

R {2) = Rl (~\ J RtBt cosv t (e iw * - 1) C[ l) + f^Bt coav t (e lwt - 1) 



[ R f B t cos v t (e 1 ^ 1 -1) -[ R tl B tl B t2 smvoh smv t 2 (l - e icJtl )ti 
Jt J Jtx.ti 



(54) 



Denoting the contribution of the last term in (54) as ^{^b) we obtain for the renormalized dissipation to 2nd order 



J_ 1 _ logi; log v-\ogv r/J_, 
r)2 V V V V2 



(55) 
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The contribution of the last term is peculiar and depends on the order of limits taken. We define a nonequilibrium 
limit where rj R is evaluated for a strictly DC field, i.e. uj — > is taken first, and then a logarithmically divergent 
E 7^ term is obtained, i.e. 



11 If 

5(-^) = — lim lim — / R tl B tl B t2 sinvh sinvt 2 (l - e lwtl )ti = 

, lim / R t .Bu sin v$ti t\ [ R t2 B t2 smvt2 = lim — =■ / sinfuti) x / sm(vt2)/tl = 

r) s v->0j tl J t2 v->0T) 6 J J " 



lim i— x v log v + 0(v) = ~ log v (56) 
v— s-o t] v rj A 

Considering next the alternative equilibrium order of limits, i.e. first E — > 0, we obtain 

lim lim sin(wti) smfv^) = (57) 

u)— fO v— ^0 



hence = 0. The renormalized 77 to second order is then 



1 1 log v log 2 v + b Q log w 



(58) 



where 60 depend on the order of limits, the nonequilibrium case has 60 = 0, while the equilibrium one has 60 = ~ 1- 
The latter case is in fact the known equilibrium result!^. The distinction between the two limits will become more 
pronounced in the full quantum treatment. 



B. Numerical solution of the Langevin Equation 

We solve the nonlinear Langevin equation numerically. The time is discretized to t = T/N x (1,2,. ..TV), with 
T the total time span of system. The noise term Q is generated numerically using a discrete Fourier transform of 
£1, = v 1 BJTW" where 1Z 1 is a unit white Gaussian noise. The correlation function linearity requires introducing a high 
frequency cutoff To- We choose the cutoff to be in Lorenzian form B u = /wy|u;|/[l + w 2 Tq], in the following section we 
explain the importance of this choice. 

We solve the equation in iterative procedure. Using the convolution form 

0t = / Rt,t> [% cos 6 t , - $ sin t , - E] (59) 



starting with an arbitrary configuration of of^ we calculate the right hand side (RHS) of (59) to find a new 0\ L> . We 
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repeat the procedure n times until the expression is saturated when 9^ — 6^ + . This procedure is improved if 
instead of taking the convolution result as the next order d t we use some mixing of that result and of the previous 
9 t configuration in the form 9^ = (1 — /3)9[ m 1 ' + j3 x RHS where /3 is mixing parameter. Typically n would be in 
order of 10 5 and /3 = 0.1. 

With this choice the Langevin equation takes the following form 

m9 t = (% cos 9 t + $ sin 9 t + E + A t (60) 
A t = 4 / sin[0 t - etAe-^^dt', 



where A t is a correction term defined by SR U 1 in the response function Eq. ( 19 ) as f t , SR t \, cos 9 t > +(^, sin 9 t > +E] = 
-J u mu 2 A u 

In the numerical system we have now four time scales, two numerical time scales, i.e. the time segment At = T/N 
and the time span T , as well as the two physical high frequency cutoffs, 1/tq for the noise and lu c the mass cutoff. The 

region of interest corresponds to velocity v R = /<? t ^ between the time scales At <C t < 1/cj c -C l/v R ~ l/v < T. 

The inequality To < 1/lu c is useful since we compare the numerical result to an asymptotic result in which u> c rather 
than 1/tq is the high frequency cutoff. 

With the result for Q t we can find the renormalized \jr\ R = dv R /dE with v R = \^t^ where the average (...) reflects 

an average on both the time domain t > 1/uj c and on numerous realizations of the noise. 
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In the left panel of Fig(2]our numerical solution for the Langevin equation is shown, including a fit to the second 
order with bo = 0. On the right panel the 1st order is subtracted with either the nonequilibrium 60 = or the 
equilibrium 6 = — 1- The first is in fact a better fit for the numerical data. When 1/v approaches the simulation time 
span T the numerics become unreliable, as the particle cannot complete even one revolution in time T; a plateau is 
then observed at low E. 




FIG. 2: Left panel: Velocity-field relation for Eq. (60 1 with 17 = 30/i/tt, oj c = 100/t and to = 20Ar. Here N = 2 1 



At = 1/20. The circles are numerical data, the full red line is a 1st order perturbation in I/77, the dashed lower red line is its 
logarithmic expansion for large lnv/ui c and the dashed upper (black) line includes the 2nd order logarithmic term, corresponding 
to Eq. |58} for b = 0. Note that the data is not reliable for E/tilu c < 1/(AtNu c ) ~ 0.06. 

E (2) _ 



Right panel: The same data and line types after subtracting the 1st order terms, i.e. 



1 



Afln JL _ 1). 



An 



additional dash-dotted line corresponds to 60 = — 1, which is a worse fit to the data than bo = (dashed upper line). Note that 
the numerical data displays E/v rather than dE/dv, hence Eq. (53 1 acquires a —1 term. 



With the numerical results for 6 T we can also generate the correlation function C T — (W T — 6*o] 2 ), the first order 
perturbation for this correlation function is given in Eq. ( |49[). In Fig. [3] we plot this correlation function as a function 
of the time separation r for the same parameters as in Fig|2J with and without a finite field. The data is fairly close 



to the 1st order result (49) for not too long times, i.e. for zero field the correlation has a subdiffusion logarithmic 
behavior while for finite force the correlation has a diffusion (~ r) behavior. 



1r 



200 r 




1000 

T 



FIG. 3: Left panel: The correlation function CV as a function of time (Blue) and the asymptotic results of Eq. (49 1 (red) for 
E = 0. Right panel: The correlation function as a function of time (Blue) and the asymptotic results of Eq. \4Q\ for Ej-q = 1 
and tq = 1. 
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IV. QUANTUM RG 
A. Perturbations from Si„t 

Consider now the definition r\ R in Eqs. fl28|29| ) 



iSE{1) = =|^5 i)t ^co S (^0sin(^ t 0cos(vi-^' + ^ t -^)^) 



2 



/ £ ; ^_ / e \iHa9 t +\iha'6 t ,+iiJ,{vt— vt'+S6 t -86 t i)\ _ 



- I Bt n a _ e -\^h{<jiR t , t -a'iR tt ,)+i^{vt-vt') /gi% 

For t < t' the term a'Rw = and then ^ a' = 0. The result is then finite only for £ > t'; defining fjf = //a', 



Hence the force correction is 



(62) 



8E {1) =~\J b t sin(|fii? T ) sin(W) (63) 



so that using Eq. (29) and performing the calculation of the integrals with arbitrary cutoffs To and ui c 1 — m/ij one 
obtains: 

4 = 1 - — [sin( A ) ln(v/u c ) + C + 0(1/*)] (64) 

where the constant C depends on tq and uj c . Although we will not need it below, its detailed form is given in the 
Appendix [C| in the limit tq = 0. 
Consider next 2nd order in Si n t, 

^ (2 M/^l\ = §4(A) 2 -A_ £ e2£3£ J Biiit2Bt3it4e ^( tl -t 2 ) +iff ',(t 3 - t4 ) 

Note that 6/S9 tl can be applied also at either £2, £3, £4 and all these terms are identical since sin(|MfJ appears in the 
same form for all ti, hence a factor 4. Now change all ej, a, a' — > — (e,, ct, ct') and define er' = er/i to obtain 



i 2 f 
18 = — ^ V] e 2 e 3 e 4 cr / B tl ,t 2 B t3jti sia[u(*i - t 2 ) + M w (*3 - 



xexp{-±h(a(E 1 § tl +e 2 d t2 +e 3 9 t3 +e 4 d t J[d tl -d t2 +i[i{9 t3 -9 t4 )}) } 



\ 12 ^£364 / B tut2 B t3iU A 2 sm{v(t 1 - t 2 + fJ.v(t 3 - t 4 )] (66) 



8h 3 



where 



A 2 = exp{i«^ei(-i? t2 , tl + M-Rta.ti _ V r um)} x 
exp{f ifte 2 (i? tl ,t 2 + M-R* 3 ,*2 ~ M-^WJ} x 
exp{|^e 3 (i? tlit3 - R t2 ,t 3 ~ nRu,t 3 )} x 

expQihei(R tl ,t 4 - Rt 2 ,t 4 + nRt 3 ,tJ} ■ ( 67 ) 
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Note that in A 2 if t 2 is the maximal time then its second factor =1 and J2 t2 £ 2 = 0- similarly, if t 3 (or t 4 ) is the 
maximal time, the the 3rd (or 4th) factor =1 and J2e 3 £ 2 = (or e 2 = 0). Therefore ii must be the maximal time 
and the 1st factor =1. The result is symmetric in t 3 ■<-» £4, so choose t 3 > t 4 , with factor 2. Hence 3 time orderings, 
denoted by A,B,C, 5E^ = SE A + 8E B + 5E C , 



SE A 



SE r 



ft 3 



E 



E 



sm[\h(R tlM + ^Rt 3 ,t 2 )} sin(|fti?t lit3 ) sin[|fi(i?t lit 



sin(f?LR tl) t 2 ) sin[|/i(i? tl]t 3 

>t 1 >t 2 >t 3 >t i 

xB tl ,t 2 Bt 3 , ti sin[u(*i - t 2 ) + fiv(t 3 - t 4 )] 
4 



*l>t3>t2>*4 

4 



tl>*3>*4>t2 



sin[|fi(i^ Iita + fiR t3it2 - (J,R ti ,t 2 )]sm(ihR tut3 ) sia[^h(R tlM + fiR t3 ,tJ] 



xB tut2 B taM sin[u(ti - t 2 ) + fJ>v(t 3 - t 4 )} 

B and C terms can be time ordered as A by t 2 <-> t 3 in B and t 2 
components, 



(68) 



^3, t 3 4-> t 2 in C. In terms of the /i = ± 



8E+ 



SEr 



8E~ A + 6EZ 



SE+ 



6E+ = 



p J siB(§ft.R tllia )sin[§fi(i2t 1) tg - Rt a ,t a )] sin[|fi(-R tl 

x[B tl , t2 B taM + B tuH B t2 ,t 3 ]sm[v(ti -h + t 3 - i 4 )] 

^3 y sin(|Si? tlit2 )sin[ifi(i?t 1)t3 - -Rt 3 ,t 3 )] sin[|S(ii tl 

x[B tu t 2 B t3 ,t 4 + B tlM B t2M ] sin[u(ti -t 2 +t 4 ~ t 3 )] 
4 

X [- B *l,t 3 - B *2,t4 



I? 



1? 



/ sin(|Si? tlit2 )sin[|S(i? tlit3 + -R t3 ,t 3 )] sin[|/i(Ji tlii4 
5 ti,i 4 - B t 2 ,i 3 ] sin[u(ti - t 3 + 1 2 - i 4 )] 



3,*4 



(69) 



In general we have two cutoffs m/rj, tq in Eq. (71 and we define T\{mjr\, To) as the cutoff time for the response Rt, 
Eq. (ItJ) . For the purpose of identifying the leading In 2 v term we take a formal limit such that this cutoff time is 
Ti — > 0. We will eventually restore physical cutoffs corresponding to m/ry, r in i? t . The only cutoff for now is r in 
B(ui), Eq. |6~|. In this limit R t — > ^Q{t)e~ St where S — > +0 to ensure the retarded nature (poles of 1/(oj + iS)). The 

0, leaving just the 

(70) 
(71) 



significant virtue of this limit is that the 1st two equations of (69) vanish since Rt lt t 3 ~ Rt 2 ,t } 
last form. The evaluation of SE^ in this limit is straightforward (Appendix D), leading to 

5E& = 4 ^sin 2 (A) sin (5) . v ln(vT )[\n(vT Q ) + 1] 



Hence from ( 29 ) 



4t: sin 2 (^) sin(-) • [ln 2 (ur ) + 31n(wr ) + 1] 



V 



So far 5E^ is calculated in a formal limit t\ — > 0. Instead of attempting the hard calculation of (69) for t\ ^ 0, 
we proceed by asserting that for any tq,t\ the leading singularity as v — > is a In 2 v term, as expected for a 2- loop 
calculation. This term must involve an r\ dependent function /^(to,ti) that has dimensions of time. Fixing the 
coefficient of ln 2 [w/ r/ (ro, t±)] as in Eq. (71 1, we have f v {jQ, 0) = tq while for tq — > 0, when t\ —tm/r]= l/u> c we must 
have the form /^(0,ri) = b{rj)T\ = b(rj)Juj c . The 2-loop correction Eq. (71) becomes at tq = 



R{2) 



4z sin 2 (f ) sin(-) • ln 2 [-6(r;)] + 0(lnv) 



The renormalized friction has therefore the form 



1 



sin — ) In — + — - sm 2 ( — ) sm(-){ln 2 — + b r? In — } 

ry 7r?7 277 w c ir z n lr\ r\ cu c cu c 



(72) 



(73) 
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We have thus identified the coefficient of the In 2 term. We expect that the In coefficient 60(77) is a smooth function of 



77 and in particular is not singular at the zeroes of sin ^j-. For h/r] <C 1 we have found 6q — 0, Eq. (58). Hence, if 60 (77) 



were singular and periodic, e.g. 60(77) ~ 1/ sin it would alter the 1st order result at h/rj <C 1, an incorrect situation. 

Furthermore a singular 60(77) is inconsistent with the RG analysis, as based on the In 2 coefficient, in subsection C. 

We note that in the semiclassical limit the perturbation expansion is in R 2n ~ 1 B n /rj 2 ~ l/?7™ +1 for large 77; in 
the quantum case the i? 2 ™ -1 factors become periodic functions. The main conclusion is that there is a new small 
parameter in the perturbation series, sin(^). 



Perturbations from S c 



Here we consider the S c interaction in Eq. ( ]21p . The S c terms are 

(§ t/ e t s c ) = (9 t ,e t s 2 c ) = o 

However, the mixed term and the corresponding correction to I/77 are 



(74) 



2 h h ft h h h ft.. 

— sin— (sin ) + — cos— (sin ) buVri) 

nn 2rj 77 77 2rj lr\ 77 77 



(75) 



which does not vanish at sin ^- = 0. Note, however that this term is ~ ft 3 , i.e. a 3 -loop term. Furthermore, other 



response functions do show such zeroes. E.g. for the Rtj 1 correlation (Eq. (77) below) we have (0t sin \k6t>S c ) = to 
1st order, while in 2nd order 



H,t' 
1 



2 . ft ft h 

— sin — Isin ) m(VTi ) 

nn rj 77 77 



(76) 



We note that there are many other operators that have vanishing perturbations at sin 

e.g. the dissipation term in Eq. ( 9 1 (9 t sm(fWt> ) ) , or the response to an AC field with frequency v (9 t cos S9f sin 1 9 



to 2nd order in Si n t, S c , 



C. RG analysis 



We note that in (731 g = ~ sin(^) acts as an unexpected small parameter for the expansion, since all divergences 



vanish when g — 0. It raises the interesting possibility that g = be viewed as a RG fixed point. For that we need 
to find a renormalized coupling which obeys multiplicative RG, the simplest choice being gn = ^ sin( 2rl R(E) )• ^ nc 
question is then whether the /3-function j3 = —EdEga. can be written only in terms of Although the non-periodic 



I/77 factor in (73) appears at first problematic, we propose that resummation from higher loops, which allows for 
higher order terms O(^) changes the 1-loop term in (73) by ~ — > sin(^). 
To further motivate this proposal we consider the response 



(77) 



Physically, e ±l 2 t> corresponds to an electric field pulse SE(t) — ±|<5(< — t') or equivalently a rapid change of flux 
by ±|, therefore Rt.t 1 correspon ds t o the difference in response to these two flux pulses. Defining the dissipation 



parameter 77^ for R t ^ as in Eq. (25) we obtain that the 1-loop term is fully periodic with 



b % . 2 h 



(78) 



hence A ^ s in(^) in Eq. p| 
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FIG. 4: RG flow of 77. 



We propose then that an RG consistent theory corresponds to 



h h 2 h , . , , 4 q, ft , , ft 

' ' ){ln'[ n wj+6o(/7jMn«J} 

(79) 



2^ = 2^ ~ n Sin2( 2^ } ln[TlW] + ^ Sin3( 2^ } C0S (^){ ln2 [ T i«] + fo o(r/) In[nt;]} 



Taking a sine of both sides it yields to order g 3 , with 60 = 60 (.9 = 0), 

9r = 9 T .9 2 ln(«/o; c ) + .g 3 [ln 2 (w/o; c ) + b a In(«/w c )] (80) 
where ± refers to g = with cos(^-) = ±1, leading to 



0Ote) = 3^ = ±5^ - b gl + 0(g%). (81) 



This RG equation is satisfied for both ± fixed points as seen by substituting (80). We propose then that g ■ = are 



exact z eroe s of the perturbation expansion and the additional requirement of an RG structure leads to the result ( 80 1 



Eq. (80) yields fixed points at J^- — nix with n = 1,2, 3, ... that are attractive at 77 > n n and repulsive at 77 < n n , 



i.e. the now of 77 7^ r\ n is always to smaller 77. At these fixed points a Gaussian evaluation yields the correlation 
(cos 9t cos$o) ~ t~ 2n . We recall now a theorem for the lattice modePSl where the equilibrium action with mass related 
cutoff is replaced by an action on a lattice resulting in an XY model with long range interactions. The theorem 
statesPl that ( cos Qt cos#q) ~ l/£ 2 ; this result was also derived^ in first order in 77. The range 77 > 771 has an RG flow 
to 771 and is therefore consistent with the theorem. The hypothesis of Gaussian fixed points corresponding to n > 2 
is inconsistent with the theorem, i.e. (cos 9 t cos 9q) becomes a relevant operator at the n < 2 points rendering them 
unstable. Note t hat in t he SEB problem cos 9 t corresponds to a lead-dot voltage and its correlations determine the 
SET conductance 1 11 ' 13 ' 20 ', while in the ring problem it corresponds to fluctuations in the circular asymmetry. 

For 77 < 771 the system could have non-gaussian fixed points or a line of fixed points as hinted by the small 77 



perturbatiorP. The equilibrium Ki((f> x ) was evaluated for small 77 and for T — >• has the form -Ki(</> x ) ~ 5{4> x — \)/T, 



i.e. the dissipation is concentrated at the single point <fi x = |. This implies from Eq. ( |39[ ) that rj R ^ T and thefore 
vanishes at temperature T = 0. It is not clear, however, that 77 = is a fixed point in the RG sense and if so what is 
its range of attraction. An 77 = fixed point would imply the implausible result that the ring conductance diverges 
for small but finite 77. We therefore expect that 771 = 77^ is the single fixed point in this system, as illustrated in Fig. 
4. 



V. DISCUSSION 



The special value 77^ = h/(2n) has a topological interpretation as a Thouless charge pump^3. Consider a slow 
change of <f> x by one unit with H<fi x — r) R (d). For the special value 77^ = H/(2tv) the total change in the position of 
the particle is J.(6)dt = 2ir, i.e. the particle comes back to the same position on the ring and a unit charge has been 
transported. Such quantization has been shown for c ases where the spectrum has a garJ^S, though quantized charge 
transport was shown also in cases without a gap 1 28 ' 29 '. The quantized 77^ also results from arguing that there should 



be a unique frequency u>E = ^-E = vasE^0 (see discussion below Eq. 38), as suggested by linear response. 
We conclude from ( 45 1 that for 77 > 771 = 77^ the SEB satisfies the quantization 



.^WiVo^A (82) 

°o e 



In particular, when rj/h> 1 we have^l from the known M*/m ~ ^l^ 1 and from Eq. Q Cq/C 9 = 1 + 0(c 7rj '/ ri ). 
We expect R q to be independent of Nq at large 77, hence 

R g = ^\l + 0(e-^ h )} (83) 
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similar to the N c = 1 case^. 

The conductance of the ring can be defined by the voltage around the ring 2irE/e and the current e(9)/2ir, hence 
we predict that the conductance for 77 > r) R is 



2 2 
e e 



Grmg ~ 4ttV ~ h • (84) 

Finally, we consider the conditions for our proposed box experiment. By analogy with E = h(f> x in the ring, we 
propose measuring the response to a gate voltage that is linear in time No — Et. This leads to a DC current into 
the Coulomb box whose dissipation is the average in Eq. ( |45"| ). The field E should be sufficiently small so that g^ is 
sufficiently near the fixed point. For an initial g ~ 1 integration of dgu/dlnE = g 2 R yields ga = 1/ \n.(hu! c /E) <C g. 
E.g. for gn < 0.1 and a typical tko c « lmeV one needs E/h < 10 8 Hz. E/h has frequency units, corresponding to 10 s 
electrons/sec flowing into the box. We propose measuring the charge fluctuations (noise) Sq(uj) — e 2 {N t N t > )u) at a 
frequency, temperature and level spacings A such that A < u, T <C 10 8 Hz, to yield the DC response of Eqs. ( 39|82 1. 
We predict then that the corresponding charge noise Sq(uj)( 2 ^) 2 A = ^ = 2tt is quantized. 
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Appendix A: Mapping the Coulomb box and the ring 



The AES mapping has been extensively used, yet we find it useful to reproduce it since the relation between 
correlation functions has received less attention. 

The Coulomb box action corresponding to the Hamiltonian ( 40 ) is 

-iHS = jl^ y d} al {ihd t - e a )d a>i - E C (N - jV ) 2 | - ihSi ead - ihS tU n 

-iHSiead = a ki( ifi dt - £fe)afe.i 

Jt k 

-iHS tun = I y^ y tk,a,ia>l i d ci + h - c - ( A1 ) 
with the partition Z = e~ s . Adding a variable 9 t to the path integral yields 

- ifiS = J\e c [N - N - ^0 t ] 2 + J2 dl&W - e a )d aji - E C {N - A ) 2 | - ihS lead - ihS tun 

= J t &dl t Sm - e a - m t )d a ,i + ^-[Mt + 2E C N } 2 ^ - ihS lead - ihS tun (A2) 
Now define d a = e~ l6t d a 

h 2 .. . „ , - 1 

ihSlead 



ihS = [IzZ dlA^dt ~ e a )d a>i + ^0 2 + 9 t N + ^ [t^al^e* 9 * + h.c] \ 

J* [ a c k,u,i J 



(A3) 



The ring action in terms of t is derived by integrating out the fermions d a and a^. Define time ordered Greens' 
functions on the dot Gn n Jul) — and on the lead Gnu = sr. In matrix notation 

GoliW) \ t ( i fc , aii e* 

"0/b,t( 



Gi%t')=[ ° Q 'o G^ t (t,f) ) + \ t* k ai e~ i(>t )*(*-0-^+T 4 (A4) 
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The trace over fermions, using det(iG) = e TrlnlG , yields 

S ef / = - E Tr ln (*> *) = - E Tr ln Ho/ (*, [*(* - *') + Goi (f, (*)] } (A5) 

Expanding in T, the Oth order is 0t independent, the 1st order vanishes, hence to 2nd order 

Seff = -\Y^Tr{G TG Q f} =-|V / G 0a ,(t,t')GokAt',t)\tk, a ,\ 2 e iet - ie t' + h.c. (A6) 

For completeness we derive the Matsubara effective action using J2 a Ga,i(T) — T n G(oj n )e lUnT with fermionic 
uj n = 7rT(2n + 1), 

G(oj n )= [ Pdot ^ = f p dot (e)[- — h - — — — ]= / Pdot(e) ? lUJr \ = -inpdot{ti)sgn(u) n ) 



Go a .i{r) = 27rp dot (0) } sin(w„r) = pdot(0) . , m \ ( A7 ) 
* — ' * — ' sm (irl t) 

where pdot{t) is the dot density of states, assumed symmetric, and eventually constant. With the lead density of 
states Piead(e); and assuming a constant f& a ^ 

S eff = -\\t\ 2 Nc Pd ot{Q)pieaM I I . 2[ ^ cos[0(r) - 6{r')] (A8) 

J J sm [7rT{T - t')\ 

where N c = is the number of channels. This is the well known equilibrium ring system with a bosonic CL 

environmentP^ where r\ = ^ir\t\ 2 N c pd o t(0)piead(0) arnd m = 1/2E C . The expansion in T is justified for \t\ 2 —> 0, 
however with N c — > oo any value of 77 can be generated. In fact N c can be fairly small and yet reproduce the N c — > 00 
case, except at exponentially small temperatures^!. A similar derivation holds for the Keldysh action leading to the 



form (21). 



We proceed now to map observables of the Coloumb box to those of the ring problem. Since the action ( A3 1 has a 



term +9 No we identify No = — <j) x where <px is the flux through the ring (in units of the quantum flux). Hence 

h{9) = f Me _ sJ' Ec (*- Ar »'5i^ +fcmio11 torms 
Je 

= [ {h6 + 2E C N - 2E c No)e~ iS ^e 2 +fcrmio„ terms = 2E C [(N) - N ] (A9) 
Je 

In particular, without interaction, tk a = 0, the charge has no fluctuations (N) = (for |jV"o| < §) so that h(9) = 
-2E c N Q =2E c <t> x . 

Consider next the time ordered T correlations (the following is the same for (dj 07,), {8^97,) with ± Keldysh 
contours), 

h 2 T{e t e\,) = [h 2 6\e t ,e-* f E ^- N °-^ 2+ic — tcrms 

Je 

= f (tW t + 2E c N t - 2E c N ){h6 t , + 2E c N t , - 2E c N )e'^^ terms 
Je 

= h 2 T{eA,)o + 4E 2 T((N t - N ){N t , - No)) (A10) 
To obtain the retarded response, 

-iV^,=9(t-t'){[A t ,B t/ }) = 9(t-t'){A t B t ,-B t ,A t ) = T{A+B+)-{B^A+) (All) 

where ± are Keldysh contour indices, so that A + is earlier than B~ . 

Define the response K t ^' of the Coulomb box, as well as the response of ring problem K t .t> in the form (displayed 
here with operators whose (A t ) = to allow relation with time ordering), 

KtM = +i9(t-t , )([(9 t -(9)),(9 t ,-(9))}) 

K tit , = +i9(t-t / )({(N t -(N)),(N t ,-(N))}) (A12) 
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Fron Eq. (AlOl we have 

h 2 T((e t - (6))(e t , - (e))) + h 2 (0) 2 = 

h 2 T(6 t 9 t >)o + 4E 2 c T((N t - (N))(N t , - (N))) + AE 2 ((N) 2 - 2N (N) + N 2 ) 



(A13) 



Now using (A9) and that the relation (AlOl holds for both terms in (All), a relation between these response functions 
is obtained 



h 2 K t +, = ~2E c h6{t-t')+4E 2 c K tit 



which is reproduced as Eq. ( 40 ) . This relation is consistent with results in Ref. 



(A14) 



Appendix B: Semiclassical case: 1st and 2nd order 
1. 1st order term 



First order perturbation of the Green's function 



Itj'l / B tlit2 {6 n 6 t2 cos(u, j (7,0 1.7, (7, 



tl,t2 



So 



—I 

T 



BtlAI ^ d ai=1 2 3 4 Exp 



iai9 tl + ia 2 t 2 + i<?9 tl - i&9 t 2 + ia 3 9 t ' + ia±9 t 



itx,n a=± 
An Averaging with Gaussian weight 

)0tl+§ t 2 + -)) 



The retarded function 



\ai=0 



K t,t> ~ 



e ivt 1 +ivt2 + ... e iRt 1 ,t 2 + iI ^t 2 .t 1 +--- 



iai(-crRt 2 ^ 1 +aiRt,t 1 )+ia2(^Rt 1 ,t 2 -oi i B.t I t 1 )+ia3(aR t t i -aR t t i+aiR, t , tl ) iav{t 1 -t 2 ) 
t\,t2 e - ' e 



hi 5>> s 

t / V] d ai B tut2 (aRt 2 .t 1 - a 4 R t ,t 1 )(^Rt 1 .t2 + )(aR tltt > - crR t2 .t' + a4-R*,ti) e 

4 J tite CT=± 

~~ / ^*i,t2 cosw(ti — t2)Rt,t 1 Rt 1 ,t 2 {Rt 1 ,t' — Rt 2 ,f) 



(Bl) 



(B2) 



iov{t\ — 1 2 ) 



In the last expression we use RtR-t = 0. 



(B3) 



2. 2nd order term 



Using the same procedure for the second order 



i 



L 1 -'tu.i 



Et u t 2 Bt 3 ,t 4 (^9t 1 0t 2 cos(9 tl — 9 t2 )9 t3 9 ti cos(9 t3 — 9 ti )9 t >9t 

*ti 7 t2,t 3 ,t 4 

x ~y ' Q ^ c ia 1 §t 1 +ia28t 2 +iot38t 3 +ia46 ti +i<7 1 9t 1 —i<Ti0t 2 +i&28t 3 —i& ^ 

CTi,o- 2 =± 



using the symmetry between <j\ f-> —o\ and t\ £2 and similarly for £3, £4 

St^B^e^* 1 -^-^-'^, \-Rt2M +Rt 3 M ~Rum +a e R tM ] 



R (2) _ 1 



' tl 5*2 5*3 )*4 

\Rti,t 2 + Rt 3 ,t 2 ~ Rt 4 ,t 2 + o^-Rt^] [-^ti,*3 — Rt 2 ,t 3 — Rt 4 ,t 3 + a6-R*,t 3 ] 
[-Rti,* 4 — Rt 2 ,t 4 + Rt 3 ,t 4 + cuaRt.ti] \Rt u t' ~ Rt2,f + Rt 3 ,t< — Rt 4 ,t' + &6Rt,t' 



(B4) 



(B5) 
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the choice t% > £2,^3,^4, only Rt.t 1 remains. R T is real, we separate the exponent to two sinus and two cosine terms 
as follow 

R f}' = \ I B tl . t2 B t3 . M (cosv(t 1 - t 2 ) cosv(t 3 - U) - sia«(*x - t 2 ) sinv{t 3 - i 4 )) i? Ml 



[Rtl,t 2 + ^t 3 ,*2 _ ^t 4 ,t 2 ] [^*1,*3 _ ^*2,*3 _ ^t4,ts] [^*1,*4 _ ^*2,«4 + Rt 3 ,t 4 ] 



(B6) 



This long multiplicity of Rt terms is now separated to 8 different terms. For the terms with the cosine we calculate 
explicitly 3 terms, which we label by a to c. Term 'a': 



B tl A 2 costal - t 2 ) x 

,*2 ,*3 i*4 

Rt,t 1 Rt 1 ,t 2 (Rt 1 ,t' — Rt 2 ,t')Bt 3 ,t 4 cosu(<3 — t4)(R tl .t 3 — Rt 2l t 3 ){Rt 1 ,t i — Rt 2 ,t 4 ) 
B tljt2 cosw(ii — t 2 )Rt,t 1 Rt 1 ,t 2 (Rt 1 ,t' — Rt 2 ,v) Ct u 



(B7) 



This term in uj space 



K = -2 R " I RtBt cosvt ^ ~ ^ Ct 



with Ct = 2(Cj1q — C\ [ ). Similarly we choose two different terms 'b' and 'c' and write them directly in 10 space 



Rl = RljRpBt cosvt (e itJt - 1) 
J R t B t cosvt (e tuJt - 1) 



R?,, — R, 



(B8) 
(B9) 



note the R^ 1 ' in the expression R b is the first order result of the retarded green function. is the reducible term 
containing multiplication of i?£ . Renormalized r\ for small v is 

h 



111 

772 2 ?7 2 
1 ft 



T}\ ITT] 2 J { 

H 2 

log V 



RtB t C{t) t 



1ix 2 r\' 

1 _ 1 

V 2 V 3 



TFT? 3 J, 

R^t Bt t = ^ 

W Jt 



, 3 logv + 0(v) 



R t B t t (log t + 7 + 0(v) + 0{l/t)) 



2n 2 rf 



log 2 v + 0(v) 



R t B t t (log t + 7 + 1 + 0(v) + 0(l/t)) 



27T 2 ?7 



h 2 r . „, m2 ft 2 



27T 2 ?7 3 



pGg« + 0(«)]' 



2n 2 rj' 



Iog 2 v + O(v) 



(BIO) 



The terms containing the sine in Eq. (B6|, are in general of order 0(v), however we have identify the following 



term which, depending on the order of limits, may contribute a term logarithmic in v for small v. 

Rt = -Rt I R tl Rt 2 B tl B t2 sin vtx smvt 2 (l - e iutl ) [ (R tl+t3 ~ R t3 ) 
Jt u t 2 Jt 3 



(Bll) 



We label the dissipation parameter form this term by S(— ^) = lim^-^oC - i^Rt an d found the logarithmic prefactor 

v 2 



in Eq. ( 56 ) , where we use for t\ > 



/ (Rt!+t 3 - Rt 3 ) = - [ (1 
Jt 3 V J-u v 



-(ti+ts)^ 



1 

V Jo 
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(B12) 
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Appendix C: Quantum case: 1st order, more details 

Let us give the detailed calculation of the first order correction in the case of a mass only cutoff, i.e. tq = 0. Taking 
the derivative of Eq. (63) in the text we have: 



i 

2 V 



TB(T)sin(-R(T))cos(vT) = — 



T>0 



— sin( — (1 — e ™ T )) cos(ut) 

r>0 T 2 V 



2-q h [ dr 



(1 — e ™ r ) cos(ut) 



— [sin(~(l - e"- T )) - sin(~)(l - e~™ T )] cos(wr)] 



r>0 



'2r? 



2-q' 



^[ sin (A)ln(-^) + /(A) + 0(l/„)] 
7r 2rj mv Lr\ 



(CI) 



since the first integral can be computed exactly and in the s eco nd one can set v = to get the constant piece. This 
determines the constant C = /(^-) given in the text in Eq. (64), where the function f(x) is defined as: 



/(*) 



+°° dt 

— [sin(x(l - e~*)) - sin(as)(l - e _t )] 

n \7~Ti r(sin(a;z) - zsinx) = -U 3 log 

o (l-z)ln(l-z) 6 



(C2) 



0(x 5 



and is a nicely convergent integral, where one can rescale t freely. Although it is not periodic in x, upon plotting it 
one notes that it seems to become almost periodic a large x. 



Appendix D: Quantum case: 2nd order for n — > 

Since sin(|/LR tl t2 ) is a retarded function, we use for R t = Q(t)e~ st 

sin[ift( J R tl)t3 + i? t2 , t3 )] -> sin(-)e- 5 ( tl - t ^- 5 (* 2 -* 3 ) 

sin[ift(J2 tllt4 - R t3 , u + Rt 2M )\ ~> sin( !L )e -Ktx-M)-6(t 3 -u)-8(t*-u) (m) 

27/ 

e.g. Fourier of t\ — t 3 and t 2 — t 3 should have l/(^i + i5)(uj 2 + iS). define the variables 

t 2 = ti — 1\, t' 3 — t 3 — t 2 , t' A = t± — t 3 

t 2 = t' 2 +t l7 h = t' 3 + t / 2 + t 1 , t A = t' i + t' 3 +t' 2 + t 1 (D2) 

These variables are more convenient since their range is independent — oo < t' 2 ,t' 3 ,t' A < 0. The product of all 
convergence factors is then e 5 ( 3t 2+ it 3+3t i ) ^ with 3,4,3 factors unimportant since 5 — > 0. Hence 

SE W = i_ sin 2 ( ^) sin A f BfyjjBfa) £ § / e-^-^) (D3) 

X r e «"l(*3+*2)+ iw 3(*4+ t 3) -|_ e i "l(*4 + *3+ t 2)+ i "2 t 3] e ,5 ( t 2+ t 3+*4) — 



— icrv + + S —iav + iwi + 5 (—2iav + iuj\ + ilj 2 + S){— icrv + iu>i + S) 

4 . 2 . H . . / ^ \ V~~* ^ / t \2 /" ^1 1 l w l| /" ^ w 2 1 l^l 

— sin ( — ) S in(-)2_, _/ (wi - <™ - z<S) 2 1 + w 2 r 2 7 ~2V w 2 - CTU - 1 + w 2 t 2 
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with the integral over u>2 
du> 2 . 



and over oj% 



vJ 2 



LU 2 



u>2 — uv — iS —U2 — crv — iS' 1 + w|tq ' ' - • w 1 - J - 

3,2 1 



00 



-av\n(vT ) + 0(v t ln(i;To) 



(D4) 



dcui I 



1 



1 



(lui — crv — iS) 2 (— lui — crv — id)' 2 1 + uj 2 Tq 
u>i 2v 2 u>i . 



duj 



)\ — v 2 (wi — crv — i5) 2 (uji + av + id) 2 



-21n(uT ) - 2 



(D5) 



where in the last integral r — > can be taken. Substituting (D4|D5) in (D3) leads to the result Eq. (70 1 
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